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This paper presents the identification of liquid slosh plant using 
the Hammerstein model based on Grey Wolf Optimizer (GWO) method. 
A remote car that carrying a container of liquid is considered as the liquid 
slosh experimental rig. In contrast to other research works, this paper 
consider a piece-wise affine function in the nonlinear function 
of the Hammerstein model, which is more generalized function. 
Moreover, a continuous-time transfer function is utilized in the Hammerstein 
model, which is more suitable to represent a real system. The GWO method 
is used to tune both coefficients in the nonlinear function and transfer 
function of the Hammerstein model such that the error between the identified 
output and the real experimental output is minimized. The effectiveness 
of the proposed framework is assessed in terms of the convergence curve 
response, output response, and the stability of the identified model through 
the bode plot and pole zero map. The results show that the GWO based 
method is able to produce a Hammerstein model that yields identified output 
response close to the real experimental slosh output. 
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1. INTRODUCTION 

Nowadays, liquid slosh inside a cargo always happens in many situations. For example, ships with 

liquid container carriers are at high risk of generating sloshing load during operation [1]. In the metal 
industries, high oscillation can spill molten metal that is dangerous to the operator [2]. Meanwhile, sloshing 
of fuel and other liquids in moving vehicles may cause instability and undesired dynamics [3]. Hence, it is 
necessary to completely study the behavior of this residual slosh induced by the container motion. One may 
study the behavior of liquid slosh through developing the exact mathematical model of liquid slosh. So far, 
many researchers focus on the first principle approach to model the slosh behavior, while there are few 
literatures to discuss it from the perspective of nonlinear system identification approach. 

On the other hand, block oriented nonlinear system identification has became a popular techniques 
to model a complex plants. The block oriented nonlinear model can be classified into three categories, which 
are Hammerstein model, Wiener model and Hammerstein Wiener model. In particular, Hammerstein model 
is a model that consists of a nonlinear function followed by linear dynamic sub-plant, while Wiener model 
consists of a linear dynamic sub-plant followed by nonlinear function, and finally, Hammerstein-Wiener 
model contains a linear dynamic sub-plant inserted between two or more nonlinear functions in series. 
Among these three block oriented models, Hammerstein model is famous due to its simple model structure 
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and it has been widely used for nonlinear system identification. Specifically, the Hammerstein model has 
been applied to model a real plant such as Solid Oxide fuel cell [4], bidirectional DC motor [5], oxygen 
uptake estimation [6], stretch reflex dynamics [7], turntable servo system [8], pneumatic muscle actuators [9], 
amplified piezoelectric actuators [10] and multi-axis piezoelectric micro positioning stages [11]. On the other 
hand, there are many tools that have been utilized to identify the Hammerstein model. There are the iterative 
method [12-14], the subspace method [15-17], the least square method [18], the blind approach [19] and 
the parametric instrumental variables method [20]. Moreover, many also consider the optimization tools for 
Hammerstein model, such as Bacterial Foraging algorithm [21], Cuckoo search algorithm [22], Particle 
Swarm optimization [23], and Genetic algorithm [24]. Based on the above literature, several limitations are 
ineluctable in their works, which are: 

a. Most of the Hammerstein models used in their study are based on discrete-time model, while many real 
plants can be easily represented in continuous-time model. 

b. Almost all the methods assume a known structure of nonlinear function, which consists of several 
basis functions. 

Though, our proposed work can solve a more general class of continuous-time Hammerstein model by 
assuming an unknown structure of nonlinear function. In particular, a piece-wise affine function is adopted 
with so many basis functions. Due to the introduction of the piece-wise affine function, a high dimensional 
design parameter tuning is considered in this study, which make the identification problem more complex. 
On the other hand, Grey Wolf Optimizer (GWO) [25] has become a top notch optimization algorithm which 
has solved various types of engineering problems [25-27]. To the best of our knowledge, there are still few 
works to discuss on the GWO for identification of Hammerstein model. Thence, it motivates us to see 
the effectiveness of the GWO in modeling the liquid slosh plant from the real experimental data. 

This paper presents the identification of liquid slosh plant using the Hammerstein model based on 
GWO method. A remote car carrying a container of liquid is considered as the liquid slosh experimental rig. 
The GWO method is used to tune both coefficients in the nonlinear function and transfer function of 
the Hammerstein model such that the error between the identified output and the real experimental output is 
minimized. The effectiveness of the proposed framework is assessed in terms of the convergence curve 
response, output response, and the stability of the identified model through the bode plot and pole zero map. 


2. RESEARCH METHOD 

This section explains on how the identification of liquid slosh based on Hammerstein model and 
GWO based method is performed. Firstly, the experimental setup of liquid slosh plant is explained. Here, it is 
shown on how the input and output data are generated. Secondly, we show the procedure to identify 
the model of liquid slosh plant using the Hammerstein based Grey Wolf Optimizer. 

2.1. Liquid slosh experimental rig 

In this study, a mobile liquid slosh plant is considered to replicate real situation of a moving 
container carrying liquid, as shown in Figure 1. In particular, a remote control car is used to carry a small 
tank filled with liquid. The tank is also equipped with four plastic wheels so that it can move smoothly as 
shown in Figure 1(a). Moreover, three accelerometer sensors (ADXL335) that are floated on the surface of 
liquid are used to measure liquid oscillation as shown in Figure 1(b). For simplicity of our study, the liquid 
slosh data from only one of the sensor is recorded and only z-axis output data is considered. Figure 2 shows 
a general schematic diagram of liquid slosh experimental rig. 



(a) 


(b) 


Figure 1. Liquid slosh experimental rig, (a) side view, (b) plan view 
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Figure 2. Schematic diagram of liquid slosh experimental rig 


In particular, an Arduino UNO is used as a data acquisition platform to process the input and output 
data. Here, we generate a voltage from the Arduino UNO to the remote car and concurrently the Arduino 
UNO also will acquire the slosh data from the accelerometer. Both the input and output data can be 
monitored and analyzed from the personal computer using the LabView software. In order to identify 
the model of liquid slosh, the remote car is required to move to a certain distance and suddenly stop to 
generate a liquid oscillation or slosh inside the tank. Thence, we apply the input voltage as shown in Figure 3 
to move the remote car. Concurrently, the liquid slosh data is recorded as shown in Figure 4. These two data 
are then used to develop the Hammerstein model based GWO, which is discussed in the next section. 




2.2. Identification of liquid slosh using Hammerstein based grey wolf optimizer 

In this section, the proposed Grey Wolf Optimizer (GWO) for identification of liquid slosh plant in 
Section 2.1 based on Hammerstein model is presented. Firstly, a problem formulation to identify the liquid 
slosh plant is explained. Then, it is shown on how to apply the GWO method to identify the liquid slosh 
based on Hammerstein model. Figure 5 shows a complete block diagram to identify the liquid slosh model in 
Section 2.1. The proposed Hammerstein model consists of nonlinear function h(u) followed by the transfer 
function G(s). The nonlinear function is a piece-wise affine function given by 



Figure 5. Block diagram of Hammerstein model based GWO 
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h(u ) = < 


c 0 + mfu - d 0 ) 
Ci + m 2 (u — d 1 ) 


vc,,-! + m a (u - dfj —!) 


d 0 <u<d 1 
d 1 < u < d 2 


if dg--! < u < do 


and the transfer function G(s) is given by 


G(s) = 


gOO 

4(s) 


s m + b m _ 1 s m - 1 4--- + b 0 

a m s m + a m _ 1 s m ~ 1 4- —I- a 0 


( 1 ) 


( 2 ) 


In (1), the symbol m ; = (cy — Cj_ 1 )/(d; — cij_ 1 )(i = 1,2, ...er) are the segment slope with connecting input 
and output points as dfi = 0,1,... cr)and q(i = 0,1,... a), respectively. For simplicity of notation, 
let d=[d o, d\, ..., do] 1 and c=[co, c\, ..., cff. The input of the real liquid slosh plant and the identified model is 
defined by u(t), while the output of the real liquid slosh plant and the identified model are denoted by y(t) 
and y(t), respectively. Thence, the expression of the identified output can be written as 

y(t) = G(s)h(u(t)) (3) 


Moreover, several assumptions are adopted in this work, which are: 

a. The order of the polynomial A(s) and B(s) are assumed to be known 

b. The nonlinear function h(u(t)) is one-to-one map to the input u(t) and the values of dfi = 1,2, ...<r) are 
pre-determined according to the response of input u(t). 

Next, let t s be a sampling time for the real experimental input and output data (u(t), y(t)) 
(t= 0, t s , 2t s , ..., Nt s ). Then, in order to accurately identify the liquid slosh model, the following objective 
function in (4) is adopted in this study: 

N 

E(G,h ) = ^(y(??t s ) -y(rjt s )) 2 (4) 

>?=o 

Finally, our problem formulation can be described as follows. 

Problem 1. Based on the given real experimental data (u(t), y(t )) in Figure 1, find the nonlinear 
function h(u) and the transfer function G(s) such that the objective function in (4) is minimized. 

Furthermore, it is shown on how to apply the GWO in solving Problem 1. For simplicity, let the 
design parameter of Problem 1 is defined as x — [b 0 b 1 ... b 7n _i a 0 % ■■■ a m c o ■■■ c a] T , 

where the elements of the design parameter are the coefficients of both the nonlinear function and the transfer 
function of the continuous-time Hammerstein model. In GWO framework, let x t (i — 1,2, ...,M) be the 
design parameter of each agent i for M total number of agents. Then, consider x t j(J = 1,2,..., D) be the /-th 
element of the vector Xj(t = 1,2, ...,M), where D is the size of the design parameter. Thence, by adopting 
objective function in (4), a minimization problem is expressed as 

“ r « £ Gw) (5) 


for iterations k= 1, 2, ..., until maximum iteration k miix . Finally, the procedure of the GWO in solving 
Problem 1 is shown below: 

Step 1: Determine the total number of agents M and the maximum iteration k m ax. Set k=0 and initialize 
the design parameter x L (0)(i = 1,2,..., M) according to the upper bound .Y up and lower bound yi ow values 
of the design parameter. 

Step 2: Calculate the objective function in (4) for each search agent i. 

- Step 3: Update the values of the best design parameter x a , the second best design parameter x /! and 
the third best design parameter x 3 based on the generated objective function in Step 2. 

- Step 4: For each agent, update the design parameter using the following equation: 


x-iC/c + 1) = 


x 1 + x 2 + x 3 
3 


( 6 ) 


where 


x 1 = x a - P ■ |Q ■ x a - Xi(k)\,x 2 = x p - P ■ |Q ■ x p - Xi(k) |, 
x 3 - x s - P ■ \Q ■ x s - Xi(k)\ 
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for z'=l, 2, M. In (6), the vectors P and Q are expressed as follows: 


P — 2a ■ r x — a 

(8) 

Q = 2 ■ r 2 

(9) 


where n and r 2 are the random vectors, where each element is generated independently from 0 to 1, and each 
element of a is linearly decreased from 2 to 0 over the course of iteration using the following equation: 


a = 



( 10 ) 


Note that the vectors of P and Q in (8) and (9), which consists of random vectors i*i and r 2 , are generated 
independently between x\ x 2 and jc 3 . 

Step 5: After the maximum iteration is achieved, record the best design parameter and obtained 
the continuous-time Hammerstein model in Figure 1. Otherwise, repeat Step 2. 


3. RESULTSAND ANALYSIS 

In this section,the effectiveness of the GWO based method for identifying the liquid slosh system 
using continuous-time Hammerstein model is demonstrated. In particular, the convergence curve response 
of the objective function in (4), the bode plot and pole-zero mapping of linear function and the plot 
of nonlinear function, will be presented and analyzed in this study. 

Based on the experimental setup in Section 2.1, the input response u(t) as shown in Figure 3 
is applied to the liquid slosh plant, and the output response y{t) is recorded as shown in Figure 4. 
Here, the input and output data are sampled at 6=0.02 for A—450. In this study, the structure of G(s) is 
selected as follows: 

B(s ) s 3 + b 2 s 2 + b 4 s + b 0 

G(s)=—r^= ---5- 5 - (11) 

A(s) a 4 s 4 + dzS 6 + a 2 s /L + ajS + a 0 

After performing several preliminary testing on the given data (u(t), j(t)). Meanwhile, the input points 
for piece-wise affine function of h(u(t)) are given by d=[ 0, 0.2, 0.4, 0.6, 0.8, 1, 2, 3, 4, 5] T . The design 
parameter x £ R 18 with its corresponding transfer function and nonlinear function is shown in Table 1. 
Next, the GWO algorithm is applied to tune the design parameter with initial values of design parameter are 
randomly selected between the upper bound x up and lower bound Xi ow as shown in Table 1. Here, we choose 
the number of agents M= 40 with maximum iterations & max = 5000. 


Table 1. Design parameter of liquid slosh plant 


X 

Coefficients 

-^low 

-^up 

x a 

X\ 

b 2 

-5 

35 

-4.8017 

x 2 

b\ 

-5 

35 

0.2657 

*3 

bo 

-5 

35 

-0.0243 

X 4 

a 4 

-5 

35 

-0.1119 

*5 

a 3 

-2200 

-1 

-4.3214 

X 5 

a 2 

-2200 

-1 

-31.6384 

X-i 

a x 

-2200 

-1 

-901.2917 

*8 

a 0 

-2200 

-1 

-1280.4483 

Xg 

Co 

-5 

5 

-1.9997 

X\o 

Cl 

-5 

5 

0.3249 

X U 

c 2 

-5 

5 

1.5968 

X \2 

C3 

-5 

5 

-0.0757 

*13 

c 4 

-5 

5 

-0.0099 

*14 

C 5 

-5 

5 

-0.0048 

*15 

C 6 

-5 

5 

-0.0251 

*16 

Cl 

-5 

5 

0.0584 

*17 

C% 

-5 

5 

0.0216 

*18 

Cg 

-5 

5 

-2.2603 
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Figure 6 shows the response of the objective function convergence after 5000 iterations to produce 
the best design parameter as shown in the final column of Table 1. It shows that the GWO based method 
is able to minimize the objective function in (4) and produce a quite close output response y(t) as compared 
to the real output y (t), which can be clearly seen in Figure 7. Note that the identified output response tends to 
yield high oscillation when input is injected to the system and it start to attenuate when the input is zero, 
which is quite similar to the response of real experimental output. 



Figure 6. Convergence curve response 


Figure 7. Response of the identified output y(t) 
and real output y(t) 


In the real experimental setup, we can say that the liquid slosh system is stable since the liquid slosh 
output is reduced gradually as t -> oo. In order to validate our model regarding the stability, we use the bode 
diagram and the pole-zero map of the identified transfer function G(s ) as shown in Figures 8 and 9, 
respectively. From the bode diagram, it can be observed that a gain margin of 2.56 dB at 17.4 rad/s and phase 
margin of 4.18 degree at 16.5 rad/s are produced, which indicates that the linear system is stable. This is also 
supported by the pole-zero map where all the poles are located at the left hand side ofy-axis. In particular, 
the obtained values of poles are -36.6623, -0.2344±/'14.5053 and -1.4828, while the obtained values 
of zeros are 4.7468, 0.0274+/0.0661. On the other hand, we also can observe the feature of nonlinear function 
by plotting the obtained piece-wise function as depicted in Figure 10. Note that our nonlinear function is not 
restricted to any form of nonlinear function (i.e., quadratic), which is more generalized and provide more 
flexibility of searching a justifiable function. 
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Figure 8. Bode plot of transfer function G(s) 
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Figure 9. Pole-zero map of transfer function G(s) 
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Figure 10. Resultant of piece-wise affine function h(u ) 


4. CONCLUSION 

In this paper, an identification of liquid slosh plant using continuous-time Hammerstein model based 
on Grey Wolf Optimizer (GWO) has been presented. The results demonstrated that the proposed generic 
Hammerstein model based on GWO has a good potential in identifying the real liquid slosh behavior. 
In particular, it is shown that the proposed method is able to produce a quite close identified output with real 
liquid slosh output. Moreover, the resultant linear model has been proved to be stable based on the bode plot 
and pole-zero map. It is also shown that the used of piecewise-affme function gives more flexibility 
for the GWO to search more generic nonlinear function. In the future, our work can be extended to various 
types of nonlinear function such as continuous-time Wiener and Hammerstein-Wiener. 
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